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Abstract 

A measurement of the underlying activity in scattering processes with a hard scale 
in the several GeV region is performed in proton-proton collisions at \/s = 0.9 
and 7 TeV, using data collected by the CMS experiment at the LHC. The produc- 
tion of charged particles with pseudorapidity \n\ < 2 and transverse momentum 
pr > 0.5GeV/c is studied in the azimuthal region transverse to that of the leading 
set of charged particles forming a track-jet. A significant growth of the average mul- 
tiplicity and scalar-pr sum of the particles in the transverse region is observed with 
increasing pj of the leading track-jet, followed by a much slower rise above a few 
GeV/c. For track-jet pj larger than a few GeV/c, the activity in the transverse region is 
approximately doubled with a centre-of-mass energy increase from 0.9 to 7 TeV. Pre- 
dictions of several QCD-inspired models as implemented in PYTHIA are compared to 
the data. 
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1 Introduction 

In hadron-hadron scatterings, the "underlying event" (UE) is defined, in the presence of a hard 
parton-parton scattering with large transverse momentum transfer, as any hadronic activity 
that is additional to what can be attributed to the hadronization of partons involved in the 
hard scatter and to related initial and final state QCD radiation. The UE activity is thus at- 
tributed to the hadronization of partonic constituents that have undergone multiple parton 
interactions (MPI), as well as to beam-beam remnants, concentrated along the beam direction. 
Good understanding of UE properties is important for precision measurements of standard 
model processes and the search for new physics at high energy. Examples are the determina- 
tion of the losses of events due to isolation criteria in lepton identification, or the computation 
of reconstruction efficiency for processes like H— > 77 where the vertex is given by the under- 
lying event. 

The first measurement of UE activity at the LHC, with proton-proton centre-of-mass energy 
a/s = 0.9 TeV, has been published by CMS [lj. The present paper, which follows the same 
analysis procedure, reports on a measurement at \/s = 7 TeV; new measurements at 0.9 TeV 
are also reported, with an event sample 30 times larger than that in Ref. [lj. In this paper all 
measurements are fully corrected for detector effects. Details are given below. The ATLAS 
collaboration has reported on measurements at v^ = 0.9 and 7 TeV [2], using slightly different 
analysis procedures. 

The UE activity at a given centre-of-mass energy is expected to increase with the hard scale in 
the interaction. Events with a harder scale are indeed expected to correspond, on average, to 
interactions with a smaller impact parameter, a feature which in turn should enhance MPI |3 . 4J. 
This increased activity is observed to reach a plateau for high scales, corresponding to an MPI 
"saturation" effect for impact parameters selected by a sufficiently hard leading interaction. 
Conversely, for events with the same hard scale but taken at different values of the centre-of- 
mass energy, MPI activity is expected to increase with y/s EHU. The present analysis is focused 
on measurements that can contribute to the understanding of the UE dynamics, through the 
comparison of events at the same y/s but with different hard scales, and the comparison of data 
with the same hard scale but with different values of y/s. 

To study the UE, it is convenient to refer to the difference in azimuthal angle, A(p, between the 
projections onto the plane perpendicular to the beam of the directions of the hard scatter and of 
any hadron in the event. With this method, the UE activity is made manifest in the "transverse" 
region with 60° < | A(p\ < 120°, even though it cannot in principle be uniquely separated from 
initial and final state radiation. In this paper, the direction of the hard scatter is identified with 
that of the leading "track-jet", i.e. the object with largest transverse momentum, pj, formed 
using a jet algorithm applied to reconstructed tracks of particles above some minimum pj 
value in the event. The leading track-jet pj is taken as defining the hard scale in the event. An 
advantage of using a track-jet as a reference is that it is an experimentally well-defined object, 
essentially free from pileup effects. No attempt is made to refer to the corresponding parton- 
level objects, as this would result in additional model uncertainties. However, the track-jet is 
much closer to the parton-level object than the leading track. Finally, in the few GeV/c region, 
the value of the track-jet pj is better defined and more stable than for calorimeter based jets, 
which suffer from large fluctuations. 

The UE dynamics are studied through the comparison with data of models implemented in 
Monte Carlo (MC) simulations adopting MPI. The predictions of the models without MPI fail 
to reproduce the evolution of the UE observables with the scale of the interaction and with 
the centre-of-mass energy 0. The predictions for inelastic events are provided here by sev- 
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eral tunes of the PYTHIA program, versions 6.420 gHH and 8.145 HZUHJ. The pre-LHC D6T 
tune |91 H0| of PYTHIA6, which describes the lower energy UA5 and Tevatron data, is a widely 
used reference that will also be used for the present analysis. The tunes DW [10J and CW [1], 
which were found to describe best the data at 0.9 TeV [1], will be discussed for the present 7 TeV 
data. The new PYTHIA6 tune, Zl |ITT| , includes pj ordering of parton showers and the new 
PYTHIA MPI model [12]. It implements the results of the Professor tunes IJT3H considering the 
fragmentation and the coulor reconnection parameters of the AMBT1 tune [14J; preliminary 
CMS UE results at 7 TeV have been used to tune the parameters governing the value and the 
\/s dependence of the transverse momentum cutoff that in PYTHIA regularizes the divergence 
of the leading order scattering amplitude as the final state parton transverse momentum pj 
approaches 0. The tune Z2 is similar to Zl, except for the transverse momentum cutoff at the 
nominal energy of y/so = 1.8 TeV which is decreased by 0.1 GeV/c. PYTHIA8 also uses the new 
PYTHIA MPI model, which is interleaved with parton showering. The new PYTHIA8 version 
4C HTBl has also been tuned to the early LHC data. The PYTHIA8 model includes soft and hard 
diffraction [16|, whereas only soft diffraction is included in PYTHIA6; the precise description 
of diffraction is, however, of little relevance for the present analyses since it has been checked 
that the diffractive contributions are strongly suppressed by the trigger and event selection 
requirements, especially for large pj values of the leading track-jet. The parton distribution 
functions (PDF) used to describe the protons are the CTEQ6L1 set [17J for D6T, Z2 and 4C, and 
CTEQ5L [18] for the other simulations. 

The outline of the paper is as follows. Section [2] presents experimental details: brief detector 
description, data samples, event and track selection, track-jet reconstruction, unfolding proce- 
dure and systematic uncertainties. Section|3]presents results on the transverse region dynamics: 
hard-scale dependence and particle spectra at a/s = 7 TeV, and centre-of-mass energy depen- 
dence of the transverse region dynamics. Section [4] summarizes the main results of the study 
and draws conclusions. 

2 Experimental Details 

A description of the CMS detector can be found in Ref. fl9|. The coordinate system has the ori- 
gin at the nominal interaction point. The z axis is parallel to the anticlockwise beam direction; 
it defines the polar angle 8 and the pseudorapidity rj = — ln(tan(0/2)). The azimuthal angle 
(p is measured in the plane transverse to the beam, from the direction pointing to the centre of 
the LHC ring toward the upward direction. The pixel and silicon strip tracker, immersed in the 
uniform 3.8 T magnetic field provided by a 6 m diameter superconducting solenoid, measures 
charged particle trajectories in the pseudorapidity range \rj \ < 2.5. The pj resolution for 1 GeV/c 
charged particles is between 0.7% at rj = and 2% at \rj\ = 2.5. 

For this analysis, the same selection conditions apply to events and tracks at 0.9 and 7 TeV; these 
conditions are very similar to those at 0.9 TeV in Ref. [1]. Minimum bias events are triggered by 
requiring activity in both Beam Scintillator Counters (BSC) Hl9ll20| , in coincidence with signals 
from both beams in the Beam Pick-up Timing for experiments (BPTX) devices |[T9"ll2"Tl ; low-p^ 
track-jets are recorded with a prescaled minimum bias trigger. At 7 TeV, in order to enhance 
the acquisition of events with a harder scale and reduce statistical fluctuations, the analysis also 
uses single-jet triggers. 

Selected events are required to contain one and only one primary vertex, reconstructed in fits 
with more than four degrees of freedom, with a z coordinate within 10 cm of the centre of the 
4 cm-wide beam collision region. Rejecting events with more than one primary vertex does not 
bias the final results, as was checked by comparing data with different pileup conditions, taken 



Table 1: Number of selected events and corresponding number of selected tracks, for three 
values of minimum track-jet pj and for both y/s = 0.9 and 7TeV. 
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leading track-jet pj > 
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3GeV/c 


20GeV/c 
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20GeV/c 


No. selected events ( x 10 3 ) 
No. selected tracks ( x 10 3 ) 


18 543 
202 952 


6 674 
120 945 


19 
638 


5140 
33 743 


783 
9 001 


0.25 
5.8 



at low and high instantaneous luminosities. 

Selected events are also required to contain a track-jet with pj > 1 GeV/c, reconstructed with 
pseudorapidity \t]\ <2. Track-jets are defined using the SISCone algorithm [22J as implemented 

7(A^ + (A 7 )2 



in the Fastjet package II23II with a clustering radius R = \/(A<p) 2 + (A?/) 2 = 0.5. Charged 
particles reconstructed in the tracker with pj > 0.5GeV/c and \rj\ < 2.5 are used to define the 
track-jet; this t] range is wider than that used for the UE analysis (\t]\ < 2) in order to avoid a 
kinematic bias. 

A track is selected for the UE analysis if it is consistent with the primary vertex and is recon- 
structed in the pixel and silicon strip tracker with transverse momentum pj > 0.5GeV/c and 
pseudorapidity \rj\ <2. A high-purity reconstruction algorithm is used, which keeps low lev- 
els of misreconstructed and poorly reconstructed tracks (24|. To decrease contamination by 
secondary tracks from decays of long-lived particles and photon conversions, the distance of 
closest approach between the track and the primary vertex is required to be less than three 
times its (significantly non-Gaussian) estimated uncertainty, both in the transverse plane and 
along the z-axis; the uncertainty in the vertex position is also taken into account. Poorly mea- 
sured tracks are removed by requiring cr(p T ) /p T < 5%, where ct(pt) is the uncertainty on the pj 
measurement. In the selected track sample with | rj \ < 2, these selections result in a background 
level of 3% : 1% from Kg and A decay products and 2% from combinatorial background, as es- 
timated using MC simulations. The number of selected events at both centre-of-mass energies, 
for track-jet pj > 1, 3, and 20 GeV/c, and the corresponding number of selected tracks are given 
in Table G] 

The distributions presented below are fully corrected for detector effects. An iterative unfold- 
ing technique [25J is used, except for some cases that will be detailed below. The PYTHIA6 
MC with Z2 tune was used to correct the experimental distributions, while Zl, D6T and the 
default configuration of pythia 8.135 ("tune 1") were used for cross-checks and systematic 
uncertainty estimates. The detector response was simulated in detail using the GEANT4 pack- 
age 1126 1 , and simulated events were processed and reconstructed in the same manner as colli- 
sion data. The simulations were found to give a very good description of all features related 
to detector performance that are relevant to this analysis. The unfolding procedure was tested 
using MC events, by comparing the genuine distributions for generated hadrons with the dis- 
tributions obtained, after unfolding, from reconstructed tracks. 

Systematic uncertainties on the corrected data have been studied in detail. They correspond 
essentially to the uncertainties described in (H taking into account the progress reported in 
Ref. [24J. They include the implementation in the simulation of vertex and track selection 
criteria, tracker alignment and tracker material content, background contamination from Kg 
and A production, trigger conditions, run-to-run variations of tracker and beam conditions, 
including the effect of pileup, and the effect of limited samples. 

Using as MC input the Zl simulation, which gives the best description of data, unfolding pro- 
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cedures were performed using both the Z2 and tune 1 models; the maximum discrepancies 
with the MC input were taken as systematic uncertainties. 

Systematic uncertainties are largely independent of one another, but they are correlated among 
data points in each experimental distribution. They are added in quadrature to statistical un- 
certainties and represented in all figures. 

3 Underlying Event in the Transverse Region 

The hadronic activity at 7 TeV in the transverse region, for charged particles with pj > 0.5 GeV/c, 
\tj\ < 2, and 60° < \A(p\ < 120°, is first presented as a function of the leading track-jet pj 
(Section |3.1| . Multiplicity and transverse momentum distributions are then reported for two 
minimal values, 3 and 20GeV/c, of the leading track-jet pj (Section [372] |. Results at the two 
centre-of-mass energies, v^ = 0.9 and 7 TeV, are finally compared (Section 
from the various pythia models are compared to the corrected data. 



3.3 '. Predictions 
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Figure 1: Fully corrected measurements of charged particles with pj > 0.5GeV/c and \t]\ < 2 
in the transverse region, 60° < \A(p\ < 120°, as a function of the pr of the leading track-jet: (left) 
average multiplicity per unit of pseudorapidity and per radian; (centre) average scalar J^pr per 
unit of pseudorapidity and per radian; (right) ratio of the average scalar 7J pj and the average 
multiplicity. Predictions of three pythia tunes are compared to the data. 

Figure [T] presents the average multiplicity and the average scalar momentum sum in the trans- 
verse region, as a function of the leading track-jet pj. For these distributions, full unfolding 
was performed in both the track-jet pj and the studied variable, for leading track-jet pj up to 
20 GeV/c. Bin-by-bin corrections were used at higher values of the leading track-jet pj where 
the pj dependence of the studied variables is small. 

The horizontal error bars indicate the bin size; the vertical inner error bars indicate the statis- 
tical uncertainties affecting the measurements; the outer error bars represent the statistical and 
systematic uncertainties added in quadrature; statistical uncertainties dominate at large values 
of the hard scale. The same conventions and considerations apply throughout this paper. 

Two components are visible for both observables in Fig. ITl a fast rise for pj < 8 GeV/c, at- 
tributed mainly to the increase of MPI activity, followed by a plateau-like region with nearly 
constant average number of selected particles and a slow increase of £pr- A similar structure is 
observed at 0.9 TeV (see Ref. [1 j and Fig. |5]below), the fast rise being limited in that case to the 
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region with leading track-jet pj < 4 GeV/c. All PYTHIA models predict such a distinct change 
of the amount of activity in the transverse region as a function of the leading track-jet pj. 

These evolutions result in a slow but continuous increase of the average pj of the selected 
particles for leading track-jet pj above a few GeV/c. This is observed in Fig. [JJ (right plot), 
which is obtained from the ratio of the two profile distributions, with the relative uncertainties 
conservatively summed in quadrature. (A similar behaviour of the ratio can be deduced at 
0.9 TeV). 
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Figure 2: Ratios, as a function of the leading track-jet pj, of three MC predictions to the fully 
corrected measurements of charged particles with pj > 0.5GeV/c and \t]\ < 2 in the trans- 
verse region, 60° < |A<^| < 120° (cf. Fig. [TJ: (left) average multiplicity; (centre) average scalar 
yj Pt) (right) ratio of the average scalar yj pj and the average multiplicity. The inner bands 
correspond to the systematic uncertainties and the outer bands to the total experimental uncer- 
tainties (statistical and systematic uncertainties added in quadrature). 

The information on the quality of the data description by the different models is summarized in 
Fig. |2J which presents the ratio of the MC predictions to the measurements for the observables 
shown in Fig. IT] Statistical fluctuations in the data induce correlated fluctuations for the various 
MC/data ratios. Variations in the error bands are related to the unfolding procedures, and to 
the different sets of data collected with different triggers. 

The description provided by Zl is very good for both the average multiplicity and the average 
scalar momentum sum, over the full leading track-jet pj range. For the PYTHIA8 4C tune, in 
the region with track-jet pj < 20GeV/c the predictions are below the data by 5% (10%) for the 
average charged multiplicity (average scalar FJpr); for larger track-jet pj values, the average 
charged multiplicity is well described but the average 7J pj is increasingly underestimated, by 
up to 20%. This confirms observations reported in Ref. ||27|. The predictions of the Z2 tune (not 
shown in this paper) reveal similar trends as for PYTHIA8 4C with, however, a more uniform 
and more limited underestimate ( < 10%) of the average yj pj. 

As illustrated by D6T, the predictions of older PYTHIA6 tunes are significantly below the data 
in the region characterized by the fast rise of the observables (track-jet pj < 8GeV/c); in the 
"saturation" region, tune D6T provides a good description of the average multiplicity but the 
average YLpr is largely underestimated for track-jet pj > 40GeV/c. In this region, DW pre- 
dictions are lower than for D6T, and CW even lower, which reflects the different values of the 
cutoff transverse momentum and its v^ dependence [1J (CW and the DW predictions are not 
shown in this paper). 

The comparison with data of the MC predictions for the ratio of the average 7J pj and the av- 
erage multiplicity is shown in Fig. [2] (right). The absolute value is described within 5% and the 
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hard-scale dependence is well described by Zl. The PYTHIA8 4C and Z2 predictions agree with 
Zl in the rising region, but the normalization is up to 10% and 20% lower in the "saturation" 
region, respectively. For D6T, the ratio is overestimated below 30 GeV/c, and underestimated 
above 30 GeV/c. 

3.2 Multiplicity and transverse momentum distributions 




-10- 4 ! 

cu : 

10" 6 C 

10 



+ Data 

— PYTHIA-6 Z1 

— PYTHIA-8 4C 

PYTHIA-6 D6T 

charged particles 

(p > 0.5 GeV/c, hi < 2, 60° < Ml < 120°) -, 

leading track-jet p > 3 GeV/c 



lL 



S10 = 



10 15 20 25 



1 


1 ' ' 


CMS\fI = 


7TeV : 


r 10" 1 






1 






*^^^ 


: 


>10" 2 






- 




: + 


Data 




^10" 3 




PYTHIA-6 Z1 
PYTHIA-8 4C 








PYTHIA-6 D6T 


- 


10" 4 


r charged particles 

= (p > 0.5 GeV/c, hi < 2, 60° < |A<f>| < 


120°) : 




" leading track-jet p > 20 GeV/c 


" 




1 i i 


l i i . i l i i i i l i i i i l i . 


i I i i i i 



10 



10 15 20 25 



O.KT 



*= 10" : 



10° 
10 7 



CMS\fs = 7TeV = 




+ Data, p T > 20 GeV/c 

+ Data, p T > 3 GeV/c 

PYTHIA Z1 , p t > 20 GeV/c 

PYTHIA Z1,p >3GeV/c 

T 

charged particles 

(p > 0.5 GeV/c, hi < 2, 60" < Ml < 120°) 



10 15 



20 



Data 

PYTHIA-6 Zl 
PYTHIA-8 4C 
PYTHIA-6 D6T 

charged particles 

(p > 0.5 GeV/c, hi < 2, 60° < Ml < 120") 

leading track-jet p > 3 GeV/c 



30 



10 B 



10 15 20 25 30 35 

Zp T [GeV/c] 



7TeV- 




10' = 



lio- 4 r 



+ Data 
— PYTHIA-6 Z1 
--- PYTHIA-8 4C 
PYTHIA-6 D6T 

charged particles 

(p > 0.5 GeV/c, hi < 2, 60° < Ml < 120") 

leading track-jet p > 20 GeV/c 



30 



10 15 20 25 30 35 

Zp T [GeV/c] 



■r 1 10 1 



Sio 2 r 



10' E 



S10 4 = 



7TeV- 




+ Data, p T > 20 GeV/c 

+ Data, p T > 3 GeV/c 

PYTHIA Z1 , p r > 20 GeV/c 

PYTHIA Z1,p >3GeV/c 

T 

charged particles 

(p > 0.5 GeV/c, hi < 2, 60" < Ml < 120°) 



25 30 



10 a 



10 15 20 25 30 35 

Zp T [GeV/c] 



£10~ i 

^10" 5 F 
-o 10 - 6 | 

10'4 
10 




I0" 1 r 

10" 2 i 



+ Data 

— PYTHIA-6 Z1 

— PYTHIA-8 4G 

PYTHIA-6 D6T 
charged particles 

(p > 0.5 GeV/c, hi < 2, 60" < M| < 120") 
leading track-jet p > 3 GeV/c 





10 12 14 

P T [GeV/c] 



10F 
1r 

Cio-'r 

ll0- 2 r 

Sio 3 r 
^io" 5 r 

T3 10 - 6 r 

io- 7 r 

10 n 




+ Data 
— PYTHIA-6 Z1 
--- PYTHIA-8 4C 

PYTHIA-6 D6T 
charged particles 

(p > 0.5 GeV/c, hi < 2, 60" < M| < 120") 
leading track-jet p > 20 GeV/c 




10 12 14 

P T [GeV/c] 



10f 

1[ 
C 10 'r 
1 10" 2 t 

o 

^10' d ^ 

"S10" ! 

"o 10" 6 k 

io-4 

10 



7TeV = 




+ Data, p T > 20 Gsyfc 

+ Data, p T > 3 GeV/c ~ 

PYTHIA Z1 , p T > 20 GeV/c 

PYTHIA Z1,p >3GeV/c 

T 

charged particles 

(p > 0.5 GeV/c, hi < 2, 60" < M| < 120") 



10 12 14 

p [GeV/c] 



Figure 3: Fully corrected measurements of charged particles with pj > 0.5GeV/c and \t]\ < 2 
in the transverse region, 60° < |A</>| < 120°: (left) normalized multiplicity distributions; (cen- 
tre) normalized scalar YIVt distributions; (right) particle pr spectra. The leading track-jet is 
required to have \tj\ < 2 and (upper row) pr > 3GeV/c, or (central row) pr > 20GeV/c. The 
plots in the lower row provide a direct comparison of the distributions for pj > 3 GeV/c and 
Pt > 20 GeV/c. Predictions of three PYTHIA tunes are compared to the data. 

Figure [3] presents, for charged particles in the transverse region, the normalized multiplicity 
distribution, the normalized 7J pj distribution, and the particle pj spectrum. Events are se- 
lected with two minimal values of the leading track-jet pr- Pt > 3GeV/c (upper row) and 
Pt > 20 GeV/c (central row). For the charged multiplicity and 7J pj distributions, full unfolding 
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was performed for leading track-jet pr > 3GeV/c, whereas for leading track-jet pr> 20GeV/c, 
in the "saturation" region, simpler unfolding is performed, which does not take into account 
the hard-scale dependence. In the latter case the unfolding procedure is found to occasionally 
introduce correlations between adjacent bins, which arise from statistical fluctuations in the 
uncorrected distributions. For the pj spectra presented in the right column of Fig. [3]bin-by-bin 
corrections were applied. The correction factors are found to be mostly independent of the 
track-jet pj and from the centre-of-mass energy. 

The distributions in Fig.|3]are presented for a range of the variables for which the total relative 
uncertainty, after unfolding, does not exceed 30%. It is remarkable that the charged particle 
spectra extend to pj > 10 GeV/c, indicating the presence of a hard component in particle pro- 
duction in the transverse region. The distributions for the two scale selections pj > 3 GeV/c and 
Pt > 20 GeV/c are directly compared in the lower-row plots of Fig-|3J Growth of the UE activity 
with increasing hard scale is observed both through multiplicity increase and single-particle pj 
spectra hardening, consistent with the increase of particle average pj shown in Fig. FT] (right). 
The three distributions are overall rather well described by the selected MC models over sev- 
eral orders of magnitude (more than 6 for the pj spectrum). 
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Figure 4: Ratios of three MC predictions to the fully corrected measurements of charged par- 
ticles with pj > 0.5GeV/c and \t]\ < 2 in the transverse region, 60° < \A(p\ < 120° (cf. Fig. [3J|: 
(left) multiplicity distributions; (centre) scalar J^pr distributions; (right) particle pr spectra. 
The leading track-jet is required to have \tj\ < 2 and (upper plots) pr > 3GeV/c, or (lower 
plots) pr > 20 GeV/c. The inner bands correspond to the systematic uncertainties and the outer 
bands to the total experimental uncertainties (statistical and systematic uncertainties added in 
quadrature). 

Detailed comparisons are provided in Fig. [4l which presents the ratio of the MC predictions to 
the measurements in Fig. [3] In the presence of a hard scale, characterized by a leading track- 
jet with p T > 20 GeV/c (lower plots in Fig. |4|, the Zl, Z2, and PYTHIA8 4C tunes describe the 
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data well in view of the steeply falling character of the distributions. They do indeed describe 
all three distributions within 10 — 15% over most of the domain, except for PYTHIA8 4C for 
very small values of N c ;, and J^pr, an d for pj > 4GeV/c. Data description by D6T is worse, 
especially the YLPt distribution and the pj spectrum. 

The description of the data in the region with leading track-jet pj > 3 GeV/c (Fig. [3]upper plots), 
dominated by interactions with a soft scale, is not so good. In this domain, all tunes overes- 
timate the contributions of events with very low multiplicity and YlPr (N c h < 4, yj pj < 
4 GeV/c); the discrepancies are largest for D6T For larger values of the observables, the predic- 
tions of Zl, Z2, and PYTHIA8 4C are reasonably close to the data, the weak points being the 
description by Zl of multiplicities between 10 and 20, and the description by all tunes of the 
Pt spectrum in the region 3 — 8 GeV/c. For D6T, as well as for DW and CW, the descriptions of 
the yj pr distribution and of the particle pj spectrum are poor. 

3.3 Centre-of-mass energy dependence 
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Figure 5: Fully corrected measurements of charged particles with pj > 0.5GeV/c and \r/\ < 2 
in the transverse region, 60° < |A<^| < 120°: (left plots) average multiplicity, and (right plots) 
average scalar yj pj, per unit of pseudorapidity and per radian, as a function of the leading 
track-jet pj, for (upper row) data at sfs = 0.9 TeV and \/s = 7 TeV; (lower row) ratio of the 
average values at 7 TeV to the average values at 0.9 TeV. Predictions of three pythia tunes are 
compared to the data. 

The centre-of-mass energy dependence of the hadronic activity in the transverse region is pre- 
sented in Fig. p| (upper plots) as a function of the leading track-jet pj, for \/s = 0.9 and 7 TeV. 
The same unfolding methodology as for Fig. Fn was applied for the data at \/s = 0.9 TeV, in 
this case with a separation between the two correction procedures at 10 GeV/c reflecting the 
narrower rising region. The large increase with \/s of the hadronic activity in the transverse 
region and its hard-scale dependence is shown in the lower plots of Fig. [5l in the form of the 



ratio of the 7 TeV to the 0.9 TeV results. Here the systematic uncertainties at 0.9 and 7TeV were 
conservatively combined quadratically thus neglecting cancellation effects. The ratios, which 
are close to one for leading track-jet pj = 1.5 GeV/c, reach a factor of two for pj > 6 — 8 GeV/c. 

The evolution with the hard scale of the ratio of the UE activity at 7 TeV and 0.9 TeV is described 
by the Zl MC. The trend is also reproduced by PYTHIA8 4C. The evolution is much too strong 
for D6T The Z2 predictions at ^/s = 0.9 TeV (not shown here) agree with Zl in shape but the 
normalization is 5-10% too high for both observables; this trend is opposite to that observed 
at 7 TeV, which indicates that a less pronounced ^s dependence of the transverse momentum 
cutoff should be adopted for tunes using the CTEQ6L1 PDF set than for tunes optimized for 
CTEQ5L. The PYTHIA8 tune 4C |T5l| already implements such a prescription. 
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Figure 6: For charged particles with pj > 0.5GeV/c and \tj\ < 2 in the transverse region, 
60° < |A0| < 120°, (left) normalized multiplicity distributions; (centre) normalized scalar TJpr 
distributions; (right) pj spectra, at \/s = 7 TeV and at \/s = 0.9 TeV. Events with leading 
track-jet pj >3GeV/c are selected. Predictions from tune Zl are compared to the data. 

The strong growth of UE activity with y/s is also striking in the comparison of the normalized 
distributions of charged particle multiplicity and of scalar YLfT as well as in the pj spectra, 
which are presented in Fig. kA for events at \/s = 7 TeV and 0.9 TeV with leading track-jet pj > 
3 GeV/c. The same unfolding methodology as for Fig. [3] (upper row) was applied at \/s = 
0.9 TeV. 



4 Summary and Conclusions 

This paper presents a study of the production of charged particles with pj > 0.5GeV/c and 
\rj\ < 2 at the LHC with the CMS detector in proton-proton collisions at \/s = 0.9 and 7 TeV. 
Events were selected according to the hard scale of the process, provided by the transverse 
momentum of the leading track-jet, which extends up to 100 GeV/c. The study was done in the 
transverse region, defined by the difference in azimuthal angle between the leading track-jet 
and charged particle directions, 60° < \Acp\ < 120°, which is appropriate for the study of the 
underlying event. All distributions were fully corrected for detector effects. 

A strong increase of the UE activity, quantified through the average multiplicity and the aver- 
age scalar transverse momentum sum of charged particles in the transverse region, is observed 
with increasing leading track-jet pj- At \/s = 7 TeV this fast rise is followed above ~ 8GeV/c 
by a "saturation" region with nearly constant multiplicity and small 7J pj increase. The large 
increase of activity in the transverse region is observed in the multiplicity distribution, in the 
yjpr distribution and in the charged particle pj spectrum, which were studied, respectively, 
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up to N c h = 30, YLVt = 35 GeV/c, and pj = 14 GeV/c. The events at the high end of the 
distributions indicate the presence of a hard component in the transverse region. 

By comparing data taken at \Js = 0.9 and 7 TeV, a strong growth with increasing centre-of-mass 
energy of the hadronic activity in the transverse region is also observed for the same value of 
the leading track-jet pj. 

The predictions of several tunes of the PYTHIA program version 6, in particular the new tunes 
Zl and Z2, and of the new version PYTHIA8 with tune 4C have been compared to the measure- 
ments. A good description of most distributions at \Js = 7 TeV and of the \Js dependence from 
0.9 to 7 TeV is provided by the Zl tune. The predictions of the Z2 and PYTHIA8 4C tunes are 
also in reasonable agreement with the data. 
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